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ABSTRACT 

About twenty years of radio observations in five bands (from 4.8 to 37 GHz) of the 
BL Lac object GC 0109+224 (S2 0109+22, RGB J0112+227), are presented and anal- 
ysed together with the optical data. Over the past ten years this blazar has exhibited 
enhanced activity. There is only weak correlation between radio and optical flares de- 
lays, usually protracted on longer timescales in the radio with respect to the optical. 
In some cases no radio flare counterpart was observed for the optical outbursts. The 
radio variability, characterised by peaks superposition, shows hints of some character- 
istic timescales (around the 3-4 years), and a fluctuation mode between the flickering 
and the shot noise. The reconstructed spectral energy distribution, poorly monitored 
at high energies, is preliminarily parameterised with a synchrotron-self Compton de- 
scription. The smooth synchrotron continuum, peaked in the near-IR-optical bands, 
strengthens the hypothesis that this source could be an intermediate blazar. More- 
over the intense flux in millimetre bands, and the optical and X-ray brightness, might 
suggest a possible detectable gamma-ray emission. 

Key words: BL Lacertae objects: general — BL Lacertae objects: individual 
(0109+224) — BL Lacertae objects: individual (0112+227) — methods: statistical 
— radiation mechanisms: nonthermal — blazars 



The compact radio-loud source GC 0109+224 (S2 0109+22, 
TXS 0109+224, RX JOl 12.0+2244, EF B0109+2228, RGB 
JOl 12+227), belonging to the Green Bank radio survey list 
C. is a synchrotron source known for more than thirty years 
JPauUny-Toth et al. Ill972ll an d was optically identified few 
years later (|Ower^^^^^on1 179771 ) . It is a strong radio- 
millimetre active galactic nucleus (AGN), with variable flux, 
polarisation degree and position angle, showing a flat aver- 
age spectrum much like the classical BL Lac objects. Unfor- 
tunately it is poorly studied beyond optical frequencies. The 
milliarcsecond (pc) scale of GC 0109+224 reveals a compact 
core, and a secondary component with no additional diffuse 
emission, and le ss luminous and /or beamed than the l-J y 
sample sources llBondi et al. II2001I: iFev fc Charlotl l200Cll . 
The kpc scale shows a faint one-sided collim ated radio jet, 
about 2 arcsec long ijWilkinson et al."lll998h . largely mis- 
aligned with the pc-scale inner region, as is frequently ob- 
served in high-luminosity and low-energy peaked BL Lac ob- 
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jects (LBLs, with the peak in infrared ba nds). This sou rce 
is a member of the 200-mJy catalog ( Marc ha et al. Il99dh . a 



sample which seems to fill the gap between the high energy 
peaked BL Lacs (HBLs, with the peak in the UV, soft-X 
bands) and the LBLs, as expected by some current blazar 
unification pictures. 

The historical optical light curve shows a behavior inter- 
mediate between the larger-amplitude variable BL Lacs, and 
the smaller-amplitude variable flat spectrum radio quasar 
(FSRQs). Increased flare activity (rapid optical variations, 
peak superposition, rapid flux drops), is clearly shown after 
1994 in the optical bands, thanks to the increased data sam- 
pling of the Perugia University Obse rvatory, which has als o 
recorded the larger and faster flares (jCiprini et al. lEo03a^ . 
During seven years GC 0109+224 showed six main opti- 
cal outbursts (weeks/month scale), and a variability mode 
placed between the flickering and the shot noise. 

Rather achromatic long term optical behaviour , is ac- 
companied by isolated outbursts with loop-like hysteresis 
of the spectral index, indicating that rapid optical variabil- 
ity is dominate d by non-thermal co oling of a single elec- 
tron population JCiprini et al. l2003al) . Both the optical flux 
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Figure 1. The complete radio-optical density flux light curves of GC 0109+224. The optical observations are historica l, and Perugia 
Obs. (Italy) data (the post-1994 data points), reported in the original or extrapolated Johnson B band (see for details ICiprini et al. I 
|22Q23)- The high sampled data at 37 and 22 GHz are from Metsahovi (Finland), and the 14.5 GHz, 8 GHz, 4.8 GHz data are from 
the long term monitoring of the UMRAO, (USA). This source clearly exhibit a higher mean-flux level, and an enhanced activity in 
the radio-mm after 1992, and in the optical after 1994, even if the optical sampling before this year is insufficient to record faster (and 
usually larger) flares. 



and the polarization are known to be variab le on differ- 
ent t i mescales. including intraday variability llsitko et al."1 
Il985l: iMead et al. 1 1199(1 : IValtaoia et aTl Il99ll) . The rela- 
tively strong variable degree (up to 30%) and direction of 
the optical linear polarisation , is one of the most noticeable 
chara cteristics of this object iXakalo |ll99lt IValtaoia et al. I 
Il993fl . even if there is no clear correlation between the 
flux level and polarisation de gree. T he host galax y of GC 
0109+ 224 is stiU unresolved (iFalomo 1996: Wright et al. I 
11998^) . but a lower lim it to the redshift z > 0.4 is sug- 
gested (|Falomolll99^ . From the old IRAS data there is 
no evidence for a thermal c omponent in the far-infrared 
lllmnev fc Neugebauer1ll988h. 

GC 0109+224 was detected in the past at X-ray bands 
by satellites Emstem, EXOSAT, ROSAT, and the source is 
a me mber of the RGB catalog ijLaurent-Muehleisen et al~l 
8. list of intermediate blazars with properties 
smoothly distributed in a large range between the LBL 
and HBL subclasses. In th e diagnostic diagram aro-cxox 
Padovani fc Giommi I Il99,'^^ GC 0109+224 appears close 
Dennett-Thorpe et al. I l200(j) to a prototype of interme- 



diate blazars like ON 231 (W Com iTosti et al. I l20o3 : 
iBottcher et al. l2002l:lTaghaferri et al. l200n|) . The HBL and 

LBL subclasses exhibit systematically distinct properties 
(i.e. Fx/Frad flux densities ratio, degree of radio core domi- 
nance, optical polarization degree and variability, etc.). The 
so called "intermediate" blazars are critical to clarify the 
relationship between these subclasses, and to validate the 
unified schemes based on bolometric power. Intermediate 
blazars often peak in optical bands (thus are selected and 
monitored mainly by optical telescopes), and during some 
variability phases the break between the two overall spec- 
tral components occurs in the X-ray bands. In some cases 
their spectral energy distributions (SEDs) can be still de- 
scribed well with synchrotron-self-Compton (SSC) models 
(like for the HBLs), and this implies, by simple scaling ar- 
guments, a possible inverse Compton (IC) emission peaked 
in GeV gamma-ray bands (see, e.g. i g.teckcr et al. 199f|). 
In the case of GC 0109+224 E GRET did not dete cted any 
gamma-ray emission (phase 1, iFichtel et al. ll]"994h . with a 
rather low upper limit of 0.01 njy above 100 MeV. 

GC 0109+224 is regularly monitored by the Univer- 
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Figure 2. The post-1994 portion of the flux curves in Fig.g] that 
match the optically best sampled period. A few observations at 
mm wavelength (87 and 250 GHz) are also added in this plot. 
In this detail radio flares clearly appear to work on considerably 
longer timescales, with respect to the optical flares. The relevant 
emission bump between the end of 1995 and mid-1999 (clearly 
visible in the 22 GHz curve), seems to be the result either of 
superposition of four flares departing from base level, or super- 
position of four peaks departing from a slower base level bump. 
Noteworthy is also the February- April 2000 outburst, because it 
seems detected (during the last decreasing phase) at all seven 
mm— radio frequencies reported in this plot. 

sity of MichiKa n Radio Astronom y Observatory (UMRAO), 
USA (see, e.g. lAller et al.1ll996l) a nd by the Metsahovi Ra - 
dio Observatory, Finland (see, e.g. iTerasranta et al. Il99jjl . 
In order to search for time correlations between the opti- 
cal and radio fluctuations at the various bands, in Section 
121 the updated radio flux data of UMRAO and Metsahovi 
are compared with t he available optical data (reported in 
ICiprini et al. ll2003a^ . In Section |^| the temporal behavior 
of the radio light curves is investigated through well known 
methods suitable for unevenly data sets. In Section |1| we 
show the reconstruction of the SED, with the available mul- 
tiwavelength data, and a first estimation of the bolometric 
energy and physical parameters using a pure SSC model. 



2 RADIO-OPTICAL CORRELATIONS 

The reconstructed optical flux history of GC 0109-1-224, and 
the last seven years of data collected mainly during the Pe- 
rugia monitoring program, are compared with the updated 
radio flux observations, belong to the Metsahovi (37 GHz 
and 22 GHz data), and UMRAO (14.5 Ghz, 8 Ghz and 4.8 
GHz data) radio observatories. At centimetre and millime- 
tre bands (up to about 90 GHz) a monthly sampling for 
blazars is in most cases enough. The source was monitored 
at an appreciable rate especially at 22 GHz, and in a signifi- 
cant long-term period at 4.8, 8 and 14.5 GHz. Some essential 
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Figure 3. Cross correlation calculat ed with the z-tra nsformed 
discrete correlation function (ZDCF, [Alexander Ill997^ . Due to 
the different timescales for the duty cycles in the radio and opti- 
cal, and to gaps empty of data in the optical, no strong evidence 
of a radio— optical correlation is found. Weak peaks in the cross- 
correlations (ZDCCF coefficient ~ 0.5 — 0.6) are found around a 
lag of 900 days (also positive in the optical-22GHz correlations, 
see text). On the other hand, the fluxes at the various radio fre- 
quencies appear well correlated, as expected, with a peak around 
the zero lag. 



properties of the centimetre e mission of this blaz ar measured 
at UMRAO are outlined in (lAller et al. I |l999h. Metsahovi 
data up to 1998 were compared with few literature optical 
observations (only 1996-1998), finding a lag of about 400 
days between optical and 22 GHz variations iHanski et al. I 
2002). 

The complete optical (Johnson B band) and radio light 
curves of GC 0109-1-224, are plotted in Fig.0 Increased ac- 
tivity was observed after 1992 in all five radio bands. Avail- 
able optical data are too few prior to 1995, to check the 
beginning of enhanced emission at this band. Since 1995 
Perugia observations have increased the sampling density, 
recording also the larger (and faster) fiares, showing well in- 
creased activity, rapid variations, and the same flares super- 
position which seems to characterise the radio flux. The flux 
history of the source, with millimetre data added, is plotted 
in Fig. 121 in more detail. Radio outbursts are clearly longer 
then the optical one, and show smaller amplitudes at lower 
frequencies. The peak frequency, in which the highest flux is 
observed during the flare, seems, at least for the outburst of 
March- Aprfl 2000 (2000.4, Fig. El, to be around millimetre 
wavelengths (as an intense flux of 3.13±0.04 and 3.31±0.06 
Jy at 250 GHz, was detec ted by the Mam bo bolometer at 
IRAM on 29 March 2000, iBertoldi lliooi ). 

The 22 and 37 GHz curves of Fig. 13 show the faster 
features in the radio, and five main peaks: around 1996.7 
(visible also at 14.5, 8 and 4.8 GHz), around 1997.5, another 
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Figure 4. The auto-correlations of the radio flux curves, calcu- 
lated with the ZDCF. These do not show any relevant feature, 
except for broad and weak {ZDACF ~ 0.6) peaks around the 
timescale of 3.14-4.06 years (at 14.5 GHz), and around 3.53-4.60 
years (at 37 GHz). From 37 GHz to 14.5 GHz, the ZDACF pro- 
file resemble to a shot noise, with a single bump in the correla- 
tion curve within the 2000 days timescale (semi-amplitude ~ 2.7 
years). 



around 1998.0 (seen also at 8 GHz), still one before 1998.5 
(seen also at 14.5 GHz, while the 37 GHz receiver was down), 
and another around 2000.4 (this is seen in all the bands). 
Other peaks are easily identifiable (Fig. around 1993.0 
(from 8 GHz to 37 GHz), and at the end of 1989 (from 
4.8 GHz to 14.5 GHz). The more violent optical activity in 
these years (Fig. |2J is characterised by six main fiares of 
much shorter duration (on the order of weeks-one month), 
with faster rise and decline. 

The noticeable temporal structure in the radio flux 
curves, between the end of 1995 and mid-1999, looks like 
an underling emission bump spreading over 3.5 years, with 
four superposed main peaks, or a superposition of four flares 
departing from the base level. It's difficult, even visually, to 
link this radio bump to a specific optical peak. The opti- 
cal emission in this period, apart of an outburst in Oct. 
1996, decreases to a ground level, with a few fiickering low 
peaks until mid-1997, then rises almost monotonically with 
a fiickering behaviour, up to beginning of 1998, and then 
displays another relevant fiare only during November 1998 
(Ciprini ct al. 2003a). Another prominent radio outburst 
occurred in February- April 2000, and seems detected at all 
the seven mm-radio frequencies (from 4.8 to 250 GHz, Fig. 
|5J|, during the final slow decreasing phase. 

Blazars usually exhibit shorter variability timescales, 
and more numerous fiares in the optical, than in mm-radio 
bands, and the optical data are often heavily affected by 
empty gaps in the observations (seasonal interruptions, and 
weather conditions). This is a problem in the search of cor- 
relations. The level of correlation can be investiga ted with 
the D i screte Correlation Function ( DCF, Edelson fc Krolik I 
Il988l: iHufnagel fc Bregmanlll992l). and the Interp olated 
Correlation Function (IGF. lOaskeU fc Peterson Il987^ ■ suit- 
able for discrete unevenly sampled data sets. In the DCF the 
number of point per time bin can vary greatly, in the IGF 
instead, the interpolation may be unreliable if the curves 
are under-sampled. We applied a method that is more ro- 



0.5 
0,0 
-0.5 
-1.0 
0.5 
0.0 
-0.5 
-1.0 

-1.5 
0.0 
-0.5 
-1.0 



-1.5 
-2.0 



r 37GHz 

R.:t,l: 


GC 0109+224 ifc, i 

,1,1111 l-i It Jjirf^''y ^ 


L 22GHz 




L 14.5GHz 




r 8GHz 

E 1 


...--}- ii-ti1-t-*1l}iP4^ ] 



2.5 3.0 
log,.(At.) [log,.(days)] 



Figure 5. The first order structure functions (SF) of the radio 
density flux light curves, in logarithmic scale. For the 37 GHz 
flux curve, the mean slope index of the steep part of the SF is 
(3 = 0.48 ±0.08, and for the 22 GHz curve f3 = 0.65 ±0.05 (bin size 
used is 40 days for both). For the 14.5 GHz series the SF mean 
slope is /3 = 0.39±0.14, and for the 8 GHz data /3 = 0.52±0.13, (50 
days the bin size used). The time scales corresponding to deepest 
drops in the SF are of 3.4-3.8 and 7.4-7.6 years. The shorter lags 
at which there are hints of a flattening in the curves, are around 
1.1-1.8 years. 

bust especially when the data are few, ensuring that there is 
a statistically meaningful number o f points in each b in: the 
Fisher z-transformed DCF, (ZDCF. lAlexander 111997^ . This 
method build data bins by equal population rather than 
equal width, and uses Montecarlo estimations for peaks and 
uncertainties. 

For a lag between 789 and 879 days (2.1-2.4 years), no 
peaks with large values were found in the optical-radio cross 
correlation curves (ZDCCF values from 0.50 to 0.62, Fig.EJ. 
Other shorter radio-optical delays for which ZDCCF peaks 
are recognizable, are around the delay of 190 days (observed 
in both the correlations with the 37 and 22 GHz fiuxes, but it 
is not meaningful because ZDCCF ~ 0.35), around a lag of 
33 days (ZDCCF ~ 0.50) between optical- 14.5 GHz curves, 
and around a lag of about 500 days {ZDCCF ~ 0.59) 
between optical-37 GHz fiuxes. Interesting to note is that 
the cross-correlations, with a peak around the 2.1-2.4 years 
lag range (both positives, delays, and negatives, leads), are 
found in all the optical-radio (37, 22, 14.5 and 8 GHz) cor- 
relation curves, (Fig. 01 . This also appears when applying 
the standard DCF method. Hints of optical-radio correla- 
tion around this lag timescale, is also present in the plots of 
iHanski et al. ll2002l) . while the 350 and 400 days lags cited 
in this work (based only on two years of optical data), was 
not confirmed by our data. 

The optical-radio delays found around 2.1-2.4 years, do 
not represent strong evidence for a real physical correlation, 
because the timescales suggested are very long. Even though 
the coefficients are not low (0.5-0.6) and the correlations are 
found for all radio frequencies correlations. On the other 
hand shorter delays around 33, 190, 500 days are not sig- 
nificant, because they are found with low value correlations 
coefficients. To investigate better the radio-optical delays, 
it is necessary to increase data sampling in future optical 
monitoring programs. 
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Figure 6. Periodogram plots of the 37, 22 and 14.5 GHz flux 
curves. Dotted lines shows the 1% false alarm significance level, 
under the hypothesis of fluctuations dominated by Poisson statis- 
tics. The timescales corresponding to the peaks above this level, 
are 1.80, 3.43, 3.74 and 4.16 years. 

3 RADIO FLUX VARIABILITY 

The radio flux auto-correlation curves (ZDACF and DACF) 
of GC 0109+224 in various bands (Fig.QJ, do not show very 
relevant feature, except for broad and weak [ZDACF ~ 0.6) 
peaks around the scale of 3.14-4.06 years in the 14.5 GHz 
curve, and around 3.53-4.60 years in the 37 GHz curve. From 
37 GHz to 14.5 GHz the auto-correlation profile resembles 
shot noise, with a single bump in the correlation curve of 
semi-amplitude ~ 2.7 years. 

In Fi e. El we plotted the first order structure fun ctions 
(SFs) re.ix.isinionetti et al. Il 98,4lH;ighes et al. Il99"^ of the 
radio light curves. The SF works in the time domain instead 
of frequencies /. It measures of the mean difference in the 
data train as a function of the separation time At in the sam- 
pling. In log-log plots the SF of an ideal time series plus mea- 
surement noise, increases and shows an intermediate steep 
curve (^^(Ai) oc (At)''). The slope P depends on the na- 
ture of the intrinsic variability, and give the Fourier power 
law index of the spectrum, (say P{f) oc 1//'^"'^^, where f is 
the frequency). The values of /3 for GC 0109-^224 (Fi^, go 
from 0.39±0.16 to 0.65±0.08, meaning an intermediate vari- 
ability mode, between the pure flickering (pink noise, /3 = 0) 
and the shot noise (red/brown noise, /3 = 1). Deep and steep 
drops in the SF plots means a little variance, and then a pos- 
sible signature of characteristic time scales: it is apparent at 
about 3.4-3.8 years and 7.4-7.6 years. Other characteristic 
timescales are given by the turnover lag where the SF profile 
begin to flatten, and in the plots of Fig. |K| these are located 
around 1.1-1.8 years, even if in our data the flattening is 
recognizable only with difficulties. 

In Fig. El are reported the Lomb-Scargle periodograms 
ijSca rglc 1982; Hornc & Baliunas 1986) of the radio curves 
which present some features. This techniq ue (calculated 
with a fast algorithm jPress fc Rvbicki Ill98gfl . analogous to 
the Fourier analysis for unevenly sampled data sets, useful 
to detect possible periodicity, and typical timescales. The 
components above the 1% false alarm threshold, annotated 
in Fig. El (from 1.8 to 4.16 years), are not very significant, 
and not easily visually identifiable in the flux light curves. 

The length of the data record, required to demonstrate 
a real recurrent timescale depends on signal-to-noise ratio, 
systematic errors, regularity in sampling, nature of the mea- 
surements, and the nature of the underlying variation. Even 



Table 1. Some roperties and fluxes of GC 0109-1-224. (1) op - 
tical counterpart Hipparcos coordinates iZacharias et al. Ill999^ : 
(2) redshift lower limit (Falomo 1996); (3) 82 c m flux den- 
sity jPouglas et al. Ill996h : (4) 20 cm flux density iOwen et al. I 
Il98nh : (5) mean flux densities from this work; (6) 87 GHz flux 
from Metsiihovi; (7) 1.2 mm Mambo-IRAM detection ('BertolchJ 
12002), (8) IRAS far-IR flux at 60 ^ira llmocy & Ncugebaue^ 
119881) : (9) mid-IR flux in L band (OdcHetal. 1978); (10) 
mean optical flux in R Cousins band iCiorini et al. 2003aJ; (11) 
mean optical flux in B Johnson band from historical light curve 
JCiprini et al. I^2003a^ : (12 ) optical (B band) m i n and max de - 
gree of liner polarization iValtaoia et al. I Il993t lTakalol lT99l: 
(13) X-ray flux at 2 ke V by Einstein Observatory (HEA0-2j IPC 
and MPC instruments <Owen et al. ll98ltlLedden fc: OdeU Il98,4 
iMaraschi et al. "1986^; (14) X-ray flux at 1 keV by CMA instru- 
mcnt on EXOSAT ( Maraschi & Maccagni1 ll988t Iciommi et al. I 
1990; Rey nolds et al. 1999); (15) X-ray flux at 1 keV b y PSPC 
instrument of RO SAT (Neumann et al. 1994; B rinkman n e^aL I 
t995; Reich et al. 2000) a nd sp ectral index i Kock et al. Ill996t 
iLaurent-Muehleisen e t al. 199 81): (16) flux lim it above 100 MeV 
by EGRET (see Tab.|5J, iFicht el et anil994) . 



Quantity 


Value 


Epoch 


Rof. 


R.A.(J2000.0)opt 


Olh 12m 05. 8238s 




(1) 


Dec.(J2000.0)opt 


+22° 44' 38.798" 




(1) 


z 


> 0.4 




(2) 


F„ (0.365 GHz) 


0.279 ± 0.021 Jy 




(3) 


F„(1.484 GHz) 


0.36 ± 0.2 Jy 


Jan 78 


(4) 


<F„(4.8 GHz)> 


0.56 Jy 


FobSO-DocOl 


(5) 


<F„(8 GHz)> 


0.67 Jy 


May78-Dec01 


(5) 


<F„(14.5 GHz)> 


0.66 Jy 


Sep79-Dec01 


(5) 


<F„(22 GHz)> 


1.07 Jy 


Jan85-DGc01 


(5) 


<F„(37 GHz)> 


1.05 Jy 


Jul84-Aug01 


(5) 


F„(87 GHz) 


2.06 ± 0.09 Jy 


2000 Apr 26 


(6) 


F^(250 GHz) 


3.31 ± 0.06 Jy 


2000 Mar 29 


(7) 


F„(60 nm) 


188 mjy 


1983 


(8) 


F^(3.5 Mm) 


11.8 mJy 


Jan77 


(9) 


<F„(0.64 /jm)> 


3.8 ± 0.2 


Nov94-Fob02 


(10) 


<F„(0.44 /jm)> 


2.4 ± 0.3 


Nov76-Fob02 


(11) 


P(S)™„/™a. 


0.7% / 29.69 % 




(12) 


F„(2 keV) 


0.15 /jJy 


Jul 1979 


(13) 


F„(l kcV) 


0.21 fiJy 


Aug 1984 


(14) 


F„(l kcV) 


0.26 ± 0.08 MJy 


Aug 1990 


(15) 


ax 


1.96 ± 0.25 


Aug 1990 


(15) 


F„(> 100 McV) 


< 0.01 njy 


1992 


(16) 



if the strength of the autocorrelations found is not high, from 
the beginning of 1990s our radio sampling is sufficient at all 
the five bands observed to consider this discovered hints of 
timescales, a starting point in the search for recurrent times 
with future observstions. As mentioned above, hints of a sort 
of periodicity (or at least of a typical timescale) around 3-4 
years, are found by applying all three methods. Future radio 
flux data, with increased sampling and extended observing 
period, and VLBI observations may or may not confirm this 
timescale. 



4 THE SPECTRAL ENERGY DISTRIBUTION 

GC 0109-1-224 is a bright X-ray point source {uF^{l keV) > 
10"^"^ erg cm~^ s~^ ), known since the observ ation of HEAO- 
1 A2 experiment (jPella Ceca et al. I [l99(1) in Jan. 1978. 
The spectral energy distribution (SED) assembled with the 
available historical multiwavelength data, is shown in Fig. 
[3 (references for data in figure's caption), and is com- 
patible with a smooth continuum given by a pure syn- 
chrotron component, like the classical BL Lac objects. In 
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Table 2. EGRET viewing periods (v. p.) and flux upper limits (u.l.) of GC 0109+224. Related quasi-simultaneous radio and optical 
fluxes from this work, are also reported in the Tabl e when available. R eferences: (a) EGRET integral flux density upper limit above 100 
MeV (first EGRET catalog, phase 1) for three v. p. iFichtel et al. Il994ft : (b) the only simultaneous optical data available iValtaoia et al. I 
Il993h . points out a low brightness and polarization degree 5.20 ±0.61; (c) u.l. given by the summed exposure s of phase 1 pointings, and 
the other EGRET v.p. of the GC 0109+224 region on phase 3 and cycle 4 (1994-95) iHartman et al. Ill999l) . Any new pointing of the 
region was done during the EGRET phase 2 (Nov. 1992-Sept. 1993, iThompson et al. Ill995^ . During EGRET cycle 8, at least one more 
pointing was available (2000, Apr. 18- Apr. 25) but our source was on the edge of the field of view. Throughout all the EGRET v.p. the 
source had a relatively low (sub-Jy) radio omission. 



EGRET 


Start-End Date 


F^{> lOOMeV) 


Ref. 


Optical Fopt 


Ref. 


Radio 




Frad 


v.p. 




[ph cm^^ s^^] 




obs. epoch [m.Jy] 




obs. epoch 




[Jy] 


0260.... 


1992 Apr 23 - Apr 28 1 










1992 Apr 26 




= 0.48 + 0.03 


0280.... 


1992 May 07 - May 14 | => 


< 6.0 X 10-** 


(a) 






1992 Apr 29 


Fir 


= 0.58 + 0.13 


0370.... 


1992 Aug 20 - Aug 27 J 






1992 Aug 28 0.94 ± 0.03 


(b) 


1992 May 06 


Fu,5 


= 0.33 + 0.01 














1992 Sep 01 


F22 


= 0.94 + 0.03 


3170.... 


1994 Feb 17 - Mar 01 ] 










1994 Mar 15 


F22 


= 0.78 + 0.11 


4250.... 


1995 Jul 25 - Aug 08 J 


< 5.6 X 10-** 


(c) 






1995 Jul 20 


F22 


= 0.75 + 0.04 














1995 Aug 04 


F22 


= 0.83 + 0.05 














1995 Aug 06 


Fir 


= 0.91 + 0.06 



Table Q we summarize properties and available multiwave- 
length fluxes of GC 0109+224. The SED data profile in Fig. 
|7| points out a synchrotron peak in the frequency range 
around the near-infred, optical (and possibly near-UV) 
range. This strengthens the assumption that GC 0109+224 
is an LBL or intermediate blazar, (as suggested for example 
in iLaur ^n^-Mi^^hlei^^n et al. Ill999l : IPennett- Thorpe et al. I 
l2000l:lBondi et al. 1120011) 

Unfortunately, a part of one rather low EGRET upper 
limit, no data is available to check the presence of the inverse 
Compton (IC) emission bump (the few and old X-ray obser- 
vations seems to have a synchrotron origin) . The 1 keV flux, 
measured by ROSAT PSPC instrument in August 1990, is 
Fx = 0.26±0.08^Jy, with ax = 1.96±0.25 (TablelHand Fig. 
0. This suggests a decreasing soft-X trend, in the vF^ rep- 
resentation. The 1 keV flux by EXOSAT CMA instrument 
in August 1984, is approximately the same {Fx = 0.21/iJy), 
while the 2keV flux by Emstem Observatory (HEAO-2) IPC 
and MPC instruments in July 1979, is Fx = 0.15 /iJy. These 
detections seem to prove the highest energy synchrotron tail 
is emitted in the soft-X ray band (as is commonly found in 
intermediate blazars). 

IRAS measured a flux density of 188 mjy at 60 iim and 
uppe r limits at the other its bands ijlmpev fc Neugebauer I 
Il988i) . As cited above GC 0109+224 seems characterized by 
a very strong (2-3 Jy) millimetre flux (a part of calibra- 
tion uncertainties). On March 29, 2000, the IRAM Mambo 
bolometer measured a flux of 3.13 and 3.31 Jy at 1.2 mm 
(250 GHz), (and aU other few available measurements are 
above 1 Jy, iBertoldi lE)oi) . While on April 26, Metsahovi 
detected a flux above 2 Jy at 87 GHz (Fi^.^and Tab.^. No 
simultaneous optical data were available during this period 
in which the source is not visible on night, but February 2000 
observations suggest a relevan t high emission state (Fi/(0.44 
/jm)~ 16 mJy). Acc ording to dValtaoia fc Terasranta Il996l : 
Ijorstad et al. ll2001^ . 7-ray flares are considered connected 
to the rising or high state of the blazar millimetre emission. 
GC 0109+224 seems on a rising phase during September 
1992 as showed by the radio flux curves, but it was not de- 
tected in 7-rays by EGRET. This could contradict the cited 
hypothesis, or the 7-rays were of short duration, or again 
probably the sensitivity was inadequate. The integral flux 
density upper limit (u.l.) above 100 MeV is rather low: the 



limit given in the flrst EGRET catalog (phase 1, Apr. 1991 
- Nov. 1992) is < 6 X 10~* photons cm~^ s~^, approxi- 
matively equivalent to 1 njy. The low optical brightness and 
the small radio fluxes (sub-Jansky), reported in Table |21 in 
correspondence with the EGRET viewing periods, point out 
a low or, in any case, non-flaring states. For details on the 
EGRET viewing periods pertaining to GC 0109+224, and 
corresponding available optical and radio fluxes, see Table 

m 

We tentatively fitted the few available data with a pure 
one-zone homogeneous leptonic SSC model, in two repre- 
sentatives (low and high) states of the source (continu- 
ous lines in Fig. |7| for the August-October 1990 quiescent 
state, and the December 2000 fiaring state). The multiwave- 
length SED data are grouped in epochs in Fig. (but we 
remark that this is only a rough incongruent and tempo- 
rally broad regroupment of the data). The applied model 
assumes a blob of dimension R embedded in a tangled 
isotropic magnetic field of mean intensity B, subjected to a 
continuous injection of shocked relativistic electrons with a 
break power law energy distribution, exponentially damped: 
Qinji'y) = Qof~^exp{-^/{k'ymax)) [cm~^ s~^], between the 
energies (Lorentz factors) 7min and "/max- This emitting re- 
gion moves in a relativistic jet, approximatively aligned with 
our line of sight and fed by an accreting supermassive black 
hole. The mechanism is described using a one-dimensional 
and time-dependent kinetic equation for the electron par- 
ticle distribution, and the ensemble synchrotron spectrum 
is convolved with the calculated distribution. The IC spec- 
trum results from the interaction of the distribution with the 
synchrotron photon field. The produced spectra were trans- 
formed to the frame of the observer (placed at an angle 6 
respect to jet direction), using the relativistic Doppler beam- 
ing bulk factor O = ((1 + z)r(l - 0cose))~\ The model is 
described for example in dCiorini fc Tostil 12003d : ICiorini I 
12003d'). The values of the parameters adopted for the two 
SED parameterisation of GC 0109+224, are reported in the 
figures caption. The steep injection (p > 2) and the exponen- 
tial damping, creates a few high energy electrons (important 
for the IC emission). Moreover the relatively high value of B 
(0.45-1.6 G), means a stronger magnetic field energy density 
Ub ~ B'^/{8tv), and a smaller radiation to magnetic energy 
ratio Urad/Us ~ Lie /Lsyn, redicing the self-Compton fiux. 
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Figure 7. The SED of GC 0109+224 assembled with the available multiwavelength data (open and filled symbols), upper limits 
(arrows), and ROSAT slope (bow-tie), at different epochs. Two SSC model fit attempts (black lines) are used, for a low emission state 
and an high state. The SED is consistent with a pure smooth synchrotron continuum, peaking between the near-IR, optical, near- 
UV, range. EGRET upper limit above 100 MeV is rather low. There is no evidence of gamma-ray emission bump, and the few X-ray 
detections suggests a downward trend. Data are grouped for no-simultaneous epochs with the same symbol. Instruments or reference s 
are in dicate in parentheses. Data from this work (UM RAO, Metsahovi and Perug ia observatories), NED database, from i ^ffiS^^^ 
198d OSC80), Mt. Lemmon jPuschell fc Stein Il98d). from JPan ll99Sl F99) and iNesci et al. l2003l N03), Einstein IPC iOwen et al 



mi. 



199C; 



1^ 



Ledden fc Odell 



1 ll98S^ . IR AS ilmoev Neugebauenil98SL IN 88'). EXOSAT CMA iMaraschi fc Maccagni Ill988t iGiommi et al, 
1999^ NOT iValtaoia et al. Ill993ll. ROSAT PSPC iNeumann et al. lll994HBrinkmann et al. lll995HKock et ah 



Reynolds et al. Ill999l. , _ _ _ . . _ 

Reich et al. Il200n^. EGRET iFichtel et al. Ill994l: iHartman et al II199I; ), 2MASS 



I ll997t) . IR AM Mambo (MarOO 



JanOl and Jan02, iBertoldi II2OO2I) and Mt. Maidanak iCiprini et al. I|2003e ). The two lines referring to the SSC fit attempts for the 
flaring state of December 2000 and the low state around the ROSAT detection on 1990, are performed with the following parameters. 
August/October 1990 low state: 'ymin = 60, "fmax = 1.7 X 10*^, "/break = 5 X 10^, pi = 2.4, p2 = 2.5, k = 0.015, B = 0.45Gauss, 



R = 9.5x lOi^'cm, 
B = 1.6Gauss, R = 



T> = 13. December 2000 high state: 7, 



8.5 X 10^'^ cm, V = 15. Redshift z = 0.4 is assumed jFalomo 11199^ 



20, 7™^=5 ><10f , -Ybreak = 7 X 10^ pi = 2.1, p2 = 2.2, k = 0.015, 



Indeed in our modelling (Fig. |7J the IC component appears 
depressed. The few multiwavelegth data available do not 
contain any hint of an high-energy spectral component. The 
values used for the physical parameters partially take in ac- 
count radio flux data in the SED. For example lower values 
of 7min could fit better the lower frequency radio data, but 
the IC emission is further depressed. The synchrotron spec- 
tral shape indicated by the available data, appears also to be 
suit able for a fit with logarithmic parabolic fun ctions (see, 
e.g. lOiomnii et al. l2nn2L-isa,mhruna et al. Probably, 
as stated for the fiux curves, the dynamics and the regions 
involved in the high energy emission are different from the 
radio ones. Anyway the usual neglect of the radio data points 
in blazar SED modelling is not ever justified. 

SED data and modelling (Fig. |7J, allude to a syn- 
chrotron peak ranging between the near-IR, for the quies- 
cent states, and the optical (perhaps also near-UV) bands 
for the flaring states. This is consistent with the view of 
GC 0109-1-224, as a blazar at the border of the LBL and 
intermediate subclasses. As a classical BL Lac, the SED is 
represented well with our pure leptonic SSC model. There- 
fore, with simple scaling considerations, the IC emission can 
be seen as the synchrotron component upshifted by jmax 
(with jmax — 10^ in this case). The peak in the opti- 
cal, in SSC descriptions, can predict a possible IC peak in 
GeV 7-ray energies through the following similarity rela- 
tion: {uGevFGev)/Lic {i'optFopt)/Lsy„, (as the the syn- 



chrotron X -ray end tail is relat ed to the TeV 7-ray tail in 
the HBLs ijStecker et al. lliAQfJ )!. Lsyn and magnetic fields 
intensity might be high for this object, reducing the IC 
dominance. However the strong millimeter flux, the peak 
of the SED in near-IR-optical bands, and the brightness 
in X-ray frequencies, might suggest the possibility of 7-ray 
emission. This emission may be detectable by the next gen- 
eration of 7-ray satellites with improved sensitivity (e.g. 
AGILE and GLAST). Moreover we note the possible cor- 
relation of GC 0109-1-224 with one of the highest energy 
cosmic-rays events ever detected. E = (1.7 — 2.6) x 10^" eV, 
on Dec. 3, 1993 jParrar fc Biermann I llOQ^) . In the sce- 
nario evading the GZK-cutoff, the primaries produced by 
particle acceleration in a blazar, should travel extragalac- 
tic distances not deflected, pointing directly to its source, 
and under this assumption the probability of a coincidental 
alignment of that event wi th GC 01094-224 was only 0.5% 
JParrar fc Biermann |[l993) . 



5 SUMMARY AND CONCLUSIONS 

In this work we have reported and analysed the largest 
amount of radio and optical data ever published on the BL 
Lac object GC 0109+224, collected over more than twenty 
years, and we have reconstructed also the most complete 
multiwavelength SED available. Some results are found, but 
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also open questions are arisen. GC 0109+224 is an exam- 
ple of a blazar which can produce strong outburst at high 
radio frequencies (2-3 Jy at > 20 GHz), but can remain 
relatively quiet (~ 0.3Jy GB6 survey) for long periods. The 
flux density variability is characterized by intermittent be- 
haviour, with a not regular alternation of relatively large 
amplitude flares, and flickering phases, both in radio and op- 
tical regimes. The source varies over all timescales sampled 
(days, months, years). Optical flares are clearly much faster 
with respect to the radio-mm bands, and some particularly 
large optical outbursts do not have obvious counterparts at 
mm and radio wavelengths. This could suggest additional 
or different emitting components responsible for the radio 
and the optical emission, and/or more rapid cooling of the 
synchrotron particles at higher frequencies. Moreover radio 
flares, due to the longer duty cycles, can blend, inhibiting 
the identification of radio counterparts of the optical flares. 

The radio light curves appear well correlated at the dif- 
ferent frequencies, around the zero lag. The radio-optical 
cross-correlation peaks, found in all the bands around the 
lag of 2.1-2.4 years, are extended and have low values (ZD- 
CCF < 0.62). The complexity, spread and low values of the 
correlation peaks, the impossibility to visually recognize the 
lags in the light curves, the very different duty-cycles, and 
especially such long delays, suggest that there is not a real 
correlation, and no physical meaning in the two-year lag. 
Other shorter delays around 33, 190, and 500 days are not 
significant because they are found with very small correla- 
tion coefficients. Fig.|21shows also that the signal is weak at 
zero time lag (between optical and 37-22 GHz curves). This 
means that strong events in the optical bands do not have 
a simultaneous recognizable counterpart in the radio bands. 
Maybe additional processes and/or difi'erent regions are re- 
sponsible for the emission at different frequencies (optical 
photons could be generated into jet components of differ- 
ent size and dynamics, in correspondence to the radio emis- 
sion). A certain level of confusion is undoubtedly created 
by induced causes (i.e. the large empty gaps in the opti- 
cal light curve data). Only with continuous monitoring over 
longer periods, with reduced gaps, will clarify the existence 
of physical radio-optical delays in this source. 

The DCF and SF shapes reflect the underlying nature of 
the process that created the radio variability. Both methods 
demonstrate, for GC 0109-1-224, a fluctuation mode between 
the flickering and the shot noise (i.e. P{f) oc 1/f"), with 
1.39 < Q < 1.65. A similar behaviour w as found also for 
the optical emission JCiprini et al. Il2003allbl) of this source. 
This power spectrum is characteristic of a random walk. In 
the case of the best sampled 22 GHz light curve, the value 
a = 1.65 found, is in strict agreement with the 5/3 value 
of the fully develop turbulence in the scalar theory of Kol- 
mogorov. Such variability could be related to fluctuations 
of the magnetic field, or of the bulk motion velocity of the 
emitting regions inside the jet. 

There is no evidence from our data, for long term peri- 
odic variations with a fixed period, but typical timescales 
in the ranges of ~ 3.2 — 4.5 years are implied, with all 
the applied methods (DCF, SF, periodogram). The more 
rapid fiickering synchrotron peaks, seem superimposed to 
long term trends in the radio. This could also be the re- 
sult of flares blending with long cooling times, as mentioned 
above. 



The increased activity and mean flux level, well ob- 
served in all the radio-optical bands, suggest that the same 
emission mechanism is responsible for radiation in both 
spectral regimes, i.e. the synchrotron radiation from shocked 
plasma into the jet, as seen also in the SED. Synchrotron 
emission in GC 01094-224 is probably dominant and pow- 
erful {Lsyn/ Lie > 1), as suggested by the high degree of 
the polarization, and by the absence of emission lines (2 
is still undetermined). The relatively high degree of lin- 
ear polarization observed could mean a weaker connection 
with the usual depolarization effects, which is a common 
affliction of blazars jets (like influence of a luminous host 
galaxy). The overall SED of GC 0109-)-224, shows a smooth 
synchrotron continuum peaked in the near-IR-optical range 
(as showed by intermediate blazars). Also our homogeneous 
SSC modelling suggests a relatively high magnetic held and 
synchrotron luminosity, diminishing the self-Compton radi- 
ation, even thought the multiwavelength data reported are 
insufficient to fully constrain the models. In particular GC 
0109+224, suffers from substantial lack of data in millimetre 
sub-mm and infrared bands, to check the importance of a 
possible thermal emission component (and also the possibil- 
ity of external-Compton contributions). 

Despite of the common deficiency of infrared data, and 
difficulties for far and mid-IR blazar monitoring, we sug- 
gest at least millimetre observations, due to the intense ra- 
diation which GC 0109+224 seems produce at these fre- 
quencies. X-ray observations of this blazar are also strongly 
encouraged, to check the presence of the high energy com- 
ponent in the spectrum. The high energy predictions for 
an intermediate blazar suffer uncertainties that become rel- 
evant in high energy tails, even when only using leptonic 
models IpBottchcr ot al. 2002). Moreover the X-ray data are 
insufficient to permit any prediction about TeV gamma-rays 
from this object (due to the suggested 2 > 0.4, TeV emis- 
sion might be interesting for studies on the extragalactic 
background light by absorption cutoffs). However millime- 
tre brightness, and synchrotron emission peaked at optical 
frequencies, could imply a GeV gamma-ray radiation de- 
tectable by the next-generation of gamma-ray space tele- 
scopes. An increased observing effort for this source, espe- 
cially at X-ray bands and beyond (for example GC 0109+224 
is just scheduled to be observed by Integral .Pian .2002 ) ) . to- 
gether with a better monitoring in radio-mm-optical bands, 
will clarify some of the questions that have arisen in our 
data and analysis. In this view our data and work will be 
useful database. 
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